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Abstract The forthcoming data expected from space mis- 
sions such as CoRoT require the capacity of the available 
tools to provide accurate models whose numerical precision 
is well above the expected observational errors. In order to 
secure that these tools meet the specifications, a team has 
been established to test and, when necessary, to improve the 
codes available in the community. The CoRoT evolution and 
seismic tool activity (ESTA) has been set up with this mis- 
sion. 

Several groups have been involved. The present paper 
describes the motivation and the organisation of this activity, 
providing the context and the basis for the presentation of the 
results that have been achieved so far. This is not a finished 
task as future even better data will continue to demand more 
precise and complete tools for asteroseismology. 
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1 Introduction 

The CoRoT satellite was launched into space on December 
27, 2006. This experiment will provide us with stellar oscil- 
lation data (frequencies, amplitudes, line widths) for stars of 
various masses and chemical compositions - mainly main 
sequence solar type stars, 5 Scuti stars and j3 Cephei stars 
- with ar i expected ac curac y on the freq uencies of a few 
10"^ Hz dBaglin et al.L l2002l : iMichel et al.ll2006h . 

Such an accuracy is needed to determine some of 
the key features of the stellar interiors such as the ex- 

tent of the convective core of intermediate mass stars 

('Roxb urgh and Voronts"ovLll999l : lMazumdar et"an.l2006l : ICunha and Metcalfel 
1 2007) or the position of the external convection zone 
JAudard and Pr ovost , 1994 ■ Monteiro. Ch ristensen-Dalsgaard and Thompsonl 
200Q : iB allot. Turck-Chieze. and GarciaLl2004j : lVerner. Chaplin and ElsworthL 



20061) and the envelope hel ium abundance in lo w mass so- 



lar ty p e stars (iMonteiro a nd Thompson, 1991; iBasu et all 
I2OO4I : IVemer. Chaplin and Elsworth. ,2006) . High-quality 
seismic data also mak e possible inverse analyses t o in- 
fer the stellar density dGough and Kosovich ev 1993 ) and 
rotati on profiles from separated multiplets (.Goupil et al.l 
1 19961) . More generally, it is expected that CoRoT will pro- 
vide constraints on those aspects of the physics of stel- 
lar interiors which ar e still poorly unde rstood, in par- 
ticular on convection dSamadietaLl l2005l) and on trans- 
port mechanisms at work in the radiative zones such as 
microscopic diffusion and rotationally induced transport 
dGoupil et all l2004l : iTheado et all l2005b . As a resuk stel- 
lar models will be improved by seismological inferences, 
and if they are combined with high-quality observational 
data of the global stellar parameters (luminosity, effec- 
tive temperature, radius, abundances), it will allow to im- 
prove the determination of those stellar parameters not di- 
rectly accessible through observation like ages and masses 
of single stars with valuable returns on the understand- 
ing of galactic structure and evolution dL ebreton, 2(3051 : 
iMonteiro. Christensen-Dalsgaard and Thompsoni . .20021) ! 

However, studies in helioseismology have taught us that 
in order to probe fine details of the solar internal struc- 
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ture, we need both the constraints of high-quality seis- 
mic data and extremely accurate numerica l solar models 
dReiter. Walsh and Weisslll995l : lGough et al.lll996l:lMonteiro. ' 



I2002I) . Therefore, to be able to draw valuable information 
from the future CoRoT data, we have to ensure that we will 
be able to interpret them with models having reached the 
optimal level of accuracy. 

In the 1990's the GONG solar model comparison project 
has been one of the first attempts to compare thoroughly so- 
lar structure models calculated with stellar evolution codes 
differing both in the numerical procedures and input physics. 
Extensive comparisons were made which allowed to better 
understand differences between mod els and to correct or im- 
prove the codes under comparison (iChr stensen-DalssaardL 
ll988HGabrielill99ll : lTurck-Chieze et al.Lll998l) . 

In the same spirit, within the CoRoT Seismology Working 
Group, the ESTA group has been set up (see Sect.O, with the 
aim to extensively test, compare and optimise the numerical 
tools used to calculate stellar models and their oscillation 
frequencies. The goals of the ESTA group have been (i) to 
be able to produce theoretical seismic predictions by means 
of different numerical codes and to understand the possible 
differences between them and (ii) to bring stellar models at 
the level of accuracy required to interpret the future CoRoT 
seismic data. 

In this introductory paper, we outline the tools and spec- 
ifications of the comparisons that will be presented in detail 
in the different papers in this volume. Sect.[2]provides a brief 
overview of the characteristics of the CoRoT seismology tar- 
gets to be modeled. Sect. [3] presents the ESTA group (par- 
ticipants, tools, past meetings and publications) and briefly 
introduces the numerical tools (stellar evolution codes and 
seismic codes) used in the successive steps of the compar- 
ison work. In Sect. IH we present the three ESTA tasks we 
have focused on. In Sect. |5] we give the specifications of 
the tasks concerning the input physics and parameters of the 
stellar models. Finally in Sect. |6] we present reference grids 
of stellar models and associated oscillation frequencies cal- 
culated by ESTA participants in parallel to the main ESTA 
work. These grids have been used all along the preparation 
of the CoRoT mission and have been made available on the 
ESTA web site. 



2 The observational seismic program of CoRoT 

The CoRoT mission has been designed and built to satisfy the 
scientific specifications as expressed in the document COR- 
SP-O-ETAN-83-PROJ, written in 2001. At that time more 
recent developments of the seismology theory as well as dis- 
coveries of pulsations were evidently not taken into account. 
More precisely CoRoT has been required to reach: 

• a photometric noise of 0.6 ppm in five days in the 0. 1 to 
10 mHz frequency domain, 

• a duty cycle of 90%, 

• a percentage of polluted frequencies less than 5% of the 
frequency interval. 



targcta_LRcO-1 
targcta_LRaQ1 
targctaJ-RcOa 
targcta_LRaQ^ 
targcta_LRc03 




log Teff. 

Fig. 1 H-R diagram of the targets of the 5 first long runs of the CoRoT 
mission, as planned presently. Long Runs are numbered as LRxOY 
where x is a or c depending on whether it is a direct run in the center or 
anticenter direction and Y is the order ot the run in the given direction. 
Standard evolutionary tracks, labelled with their mass (in solar units) 
and corresponding to solar metall icity, no overshooting, are also shown 
(see lLebreton and Michelll2008al) . 



The major goals of the seismology programme were then 
to detect solar type oscillations in a few bright objects, cov- 
ering the H-R diagram and to see whether classical pulsators 
(5 Scuti and j3 Cephei stars in particular) near the main se- 
quence had also lower amplitude modes. With these goals 
in mind the seismology program has been defined with two 
components: 

★ the central program having at least 5 sessions of 150 days 
(see Fig.[T|, 

★ the exploratory program with a few short sessions of 20 
days, to complement the H-R diagram coverage. 

The instrument has two different fields: 

• the seismology field in which 10 bright stars are ob- 
served at the same time (5.5<mv<9.5) with a sampling 
of 32 seconds, 

• the exoplanet field which observes 12 000 stars at the 
same time (11 Kiny < 16) with a sampling of 5 12 seconds 
and in a few cases 32 seconds. In this field the required 
photometric accuracy is 7x10"^ in 1 hour for a 15.5 
magnitude star. 

After one year of operations the mission fulfils all its 
specifications. 



Goals and Tasks of ESTA-CoRoT 



3 



3 ESTA organisation and tools 

The Evolution and Seismic Tools Activity (ESTA) has been 
set up in 2002, following the ES A announcement of opportu- 
nity for participation in the CoRoT mission. The actual work 
has started in 2004 at the kick-off meeting that took place 
during the CoRoT Week 7 in Granada (December 2004). 
Following this several other meetings and workshops (7 so 
far) have been organised as the activities progressed. A web- 
site has been created to coordinate the activities, allowing 
for an efficient exchange of data and documentation as well 
as to share the results of the tasks that have been organised 
around the key objectives of ESTA. This site is available at 
the following url: 

http : //www. astro. up.pt/corot/ 

ESTA has been open to participation of the whole com- 
munity and in particular the involvement of groups with evo- 
lution codes or seismic codes has been actively pursued. Up 
to now there are over 40 team members involved in one 
or several of the tasks corresponding to more than 16 in- 
stitutes/groups spread over 9 countries. It includes the par- 
ticipation of 10 evolution codes and 9 oscillation codes in 
the tasks and comparisons that have been developed over 
the last three years. As part of the effort several results have 
been published (as papers or proceedings of the workshops) 
or made available through the web-page. 



3.1 Stellar internal structure and evolution codes 

Ten stellar evolution codes have participated in ESTA at dif- 
ferent levels. The list is given below. Each of them is pre- 
sented in detail in a specific paper in this volume. 

• ASTEC - The Aarhus Stellar Evolution Code is a Danish 
stellar evolution code that has b een involved in all the 
ESTA T ASKS. It is presented by IChristensen-DalsgaardI 
(I2008bh . 

• ATON - This is an Italian stellar evolution code. Up to 
now it has participated onl y marginally in ESTA activi- 
ties. It is presented by Ventu ra et alj ((2008). 

• CESAM - The Code d' Evolution Stellaire Adaptatif et 
Modulaire is a French stellar evolution code. This public 
cod^] has been involved in all the ESTA TASKs as well 
as in the reference grids c alculation. It is presented by 
iMorel and LebretonI (l2008h . 

• CLES - The Code Liegeois d' Evolution Stellaire is a Bel- 
gian stellar evolution code that has been involved in all 
the ESTA TASKs as well as in the reference grids calcu- 
lation. It is presented bv lScuflaire et a l. (2008b). 

• FRANEC - The Frascati Raphson Newton Evolutionary 
Code is an Italian stellar evolution code that has par- 
ticipated in ES TA TASKs 1 and 3. It is presented by 
iDeglTnnocentil &lOi) . 



• GARSTEC- The Garching Stellar Evolution Code is a 
German stellar evolution code that has participated in 
ESTA T ASKS 1 and 3. It is presented bv lWeiss and Schlattll 
(I2008h . 

• GENEC - The Geneva Stellar Evolution Code is a Swiss 
stellar evolution code t hat has participate d in E STA 
TASK 1. It is presented bv lEggenberger et all ( |2008|) . 

• STAROX - This is an English stellar evolution code that 
has partic ipated in ESTA TASK 1. It is presented by 
iRoxburghl (io08b). 

• TGEC - The Toulouse-Geneva Evolution Code is a 
French stellar evolution code that has participated in 
ESTA TASKS 1 and 3. It is presented by iHui-Bon-Hoal 
(l2008h . 

• YREC - The Yale Rotating Stellar Evolution Code is a 
US stellar evolution code. Up to now it has participated 
only marginally in ES TA activities. It is presented by 
iDemarque et alJ (l2008h . 



3.2 Stellar oscillation codes 

Nine stellar oscillation codes have been used for the compar- 
isons of ESTA TASK 2. The list is given below. Each of them 
is presented in details in a specific paper in this volume. 

• ADIPLS - The Aarhus Adiabatic Pulsation Packag^ 
contains the ADIPLS stellar oscillation code itself but also 
many programs to manipulate the stellar model files and 
the output o f the pulsation program (cf. S ect. |5]). It is 
presented by IChristensen-Da lsgaardl (l2008ah . 

• FILOU - This stellar oscillations code developed at Paris- 
Meudon Observatory is presented bv lSuarez and Goupill 
(I2008h . 

• GRACO - The Granada O scillation Code is presented by 
iMova and Garridol (l2008h . 

• LOSC - The Liege O scillations Code is presented by 
IScuflaireetal.l(l2008ah . 

• NOSC - The Nice Oscillations Code is presented by 
iProvostI ( 120081) . 

• OSCROX - This English stellar oscillations code is pre- 
sented by Roxburgh (2008a). 

• POSC - The Po rto Oscillations Code is presented by 
iMonteirol (|2008|) . 

• PULSE - This Canadian stellar oscillatio ns code is pre- 
sented by lBrassard and Charpinetl (l2008h . 

• LNAWENR - This Romanian code for Linear Non- 
Adiab atic NonRadial Waves is presented by ISuranI 
(I2008h . 



4 Presentation of tlie TASKs 

To perform the comparison of the participating numerical 
codes, the ESTA group has focused on specific TASKs. The 



available at |http : //www . obs-nice . f r/ cesam] 



^ available at |http: //astro. phys . au. dk/~jcd/adipack.n| 
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first task, TASK 1, has consisted in comparing stellar mod- 
els and evolution sequences produced by eight stellar in- 
ternal structure and evolution codes. For that purpose we 
have fixed some standard input physics and initial param- 
eters for the stellar models to be calculated and we have 
defined several specific study cases corresponding to stars 
covering the range of masses, evolution stages and chemi- 
cal compositions expected for the bulk of CoRoT target stars. 
The second task, TASK 2, has consisted in testing, compar- 
ing and optimising the seismic codes by means of the com- 
parison of the frequencies produced by nine different oscil- 
lation codes, again for specific stellar cases. Finally, the third 
task, TASK 3, still in progress, is similar to TASK 1 but for 
stellar models including microscopic diffusion of chemical 
elements. 



4.1 TASK 1: basic stellar models 

TASK 1 has been defined with the aim of comparing the 
properties of stellar models and evolutionary sequences cov- 
ering a rather large range of stellar parameters (initial mass, 
initial chemical composition and evolution stage) corre- 
sponding to the bulk of the CoRoT targets. 

Seven study cases have been considered, the specifica- 
tions of which are given in Table [T] The cases are defined 
as follows. Seven evolutionary sequences have been cal- 
culated for different values of the stellar mass and initial 
chemical composition (X,Y,Z where X, Y and Z are re- 
spectively the initial hydrogen, helium and metallicity in 
mass fraction). The masses are in the range 0.9— S.OMq. 
For the initial chemical composition, several pairs {Y,Z) 
have been considered by combining two different val- 
ues of Z (0.01 and 0.02) and two values of Y (0.26 
and 0.28). The corresponding values of (Z/X) are in the 
range 0.014-0.029 We have adopted the so-c alled GN93 
solar mixture (see iGrevesse and NoelsL Il993h which has 
{Z/X)c.) = 0.0245. We have therefore metallicities such that 
[Fe/H] e [-0.25, +0.17] dex (with the standard definition 
[Fe/H]= log (Z/X)- log (Z/X)q). 



Table 1 Target models for TASK 1 (see Fig.[2j. We have considered 7 
cases corresponding to different initial masses, chemical compositions 
and evolutionary stages. One evolutionary sequence (denoted by "Ov" 
in the 5* column has been calculated with core overshooting - see 
text). 



Case 


M/Mv: 


Yo 


Zo 


Specification 


Type 


1.1 


0.9 


0.28 


0.02 


Xc=0.35 


MS 


1.2 


1.2 


0.28 


0.02 


Xc=0.69 


ZAMS 


1.3 


1.2 


0.26 


0.01 


Mf^<==O.lOM0 


SGB 


1.4 


2.0 


0.28 


0.02 


rc = 1.9xl0^ K 


PMS 


1.5 


2.0 


0.26 


0.02 


Xc=0.01,Ov 


TAMS 


1.6 


3.0 


0.28 


0.01 


Xc=0.69 


ZAMS 


1.7 


5.0 


0.28 


0.02 


Xc=0.35 


MS 




Fig. 2 H-R diagram showing the targets for TASK 1 (see Table [TJ. 
Red lines correspond to the PMS, black lines to the MS and blue lines 
to the post-main sequence SGB evolution. The targets are ordered in 
mass and age along the diagram, from Case 1 . 1 (bottom-right) up to 
Case 1.7 (top-left) 



On each evolutionary sequence, an evolutionary stage 
has been selected, either on the pre-main sequence (PMS), 
main sequence (MS) or subgiant branch (SGB). On the PMS 
we have specified the value of the central temperature of the 
model (rc = 1 .9 X lO'' K). On the MS, we have fixed the value 
of the central hydrogen content: X^ = 0.69 for the model 
close to the zero age main sequence (ZAMS), X^. = 0.35 for 
the model in the middle of the MS and X^ = 0.01 for the 
model close to the terminal age main sequence (TAMS). On 
the SGB, a model can be chosen by specifying the value 
of the mass M^^ of the central region of the star where 
the hydrogen abundance is such that X < 0.01. We chose 
M^" = 0.10 Mq. Figure |2] displays the location in the H-R 
diagram of the targets for TASK 1 together with their parent 
evolutionary track. 

All models calculated for TASK 1 are based on rather 
simple input physics, currently implemented in stellar evo- 
lution codes and one model has been calculated with over- 
shooting (see Sect.|5]below). 

The first results of TASK 1 have been presented by 
iMonteiro et al.' (2006). Detailed results are p resented in this 
volume (Cebreton et al.. .2008h : .Montalbdn et al.L l2008bl: 
iMarconi et al.ll2008h . 



Goals and Tasks of ESTA-CoRoT 



5 



4.2 TASK 2: oscillation frequencies 

Under TASK 2 we have aimed at evaluating the precision 
and uncertainties affecting the determination of oscillation 
frequencies. In particular we have studied the numerical pre- 
cision of the frequencies using different meshes and with a 
wide range of options to solve the equations for linear adi- 
abatic (radial and non-radial) oscillations. In particular two 
steps have been considered in this TASK: 

Step 1 - comparison of the frequencies from different 
codes for the same stellar model; 

Step 2 - comparison of the frequencies from the same 
seismic code for different stellar models of the same stel- 
lar case. 

In this exercise stellar models at the middle of the main se- 
quence, produced by ASTEC and CESAM, have been used 
with a mass of 1 .5 Mq. The models have been provided with 
different number of mesh points. 

Pr eliminary results of this task were reported by iMoval 
(1200 7') while the final result s of Step 1 are discussed in this 
volume by .Moya et al.l (l2008i ). Step 2 remains to be done. 



4.3 TASK 3: stellar models including microscopic diffusion 

The goals of TASK 3 are to test, compare and optimise 
stellar evolution codes which include microscopic diffusion 
of chemical elements. At this stage we only consider dif- 
fusion resulting fr om pressure, temperature an d concentra- 
tion gradients (see lThoul and Montalbanll2007l) while we do 
not take into account diffusion due to the radiative forces, 
nor the extra-mixing of chemical elements due to differen- 
tial rotation or inte r nal gr a vity w aves (see lAlecianL l2007l : 
iMathis et all l2007l : IZahnl |2007|) . The other physical as- 
sumptions proposed as the reference for the comparisons are 
the same as used for TASK 1 and no overshooting (see Sect. 

m. 

Three study cases have been considered for the mod- 
els to be compared. Each case corresponds to a given value 
of the stellar mass (see Table |2|. We chose rather low val- 
ues of the masses (i.e. M < 1.4 Mq) in order to keep in 
a mass range where radiative accelerations can be omitted. 
Furthermore, this avoids the problems occurring at higher 
masses where the use of microscopic diffusion alone pro- 
duces a very substantial depletion of helium and heavy ele- 
ments at the surface (and a concomitant increase of the hy- 
drogen content) and in turn requires to invoke other mixing 
proces ses to control the gr avitational settling (see, for in- 
stance iTurcotte et all Il998h . 

For the three cases we have adopted a chemical compo- 
sition close to the solar one {Z/X = 0.0243). For each case, 
models at different evolutionary stages have been consid- 
ered. We focused on three particular evolution stages: mid- 
dle of the MS, TAMS and SGB (respectively stage A, B and 
C). Figure [3]displays the location in the H-R diagram of the 
targets for TASK 3 together with their parent evolutionary 
track. 



Table 2 Target models for TASK 3. Left: Three cases with corre- 
sponding masses and initial chemical composition. Right: Three evolu- 
tionary stages examined for each case. Stages A and B are respectively 
in the middle and end of the MS stage. Stage C is on the SGB. 



Case 




Yo 


Zo 


Stage 


X, 




3.1 


1.0 


0.27 


0.017 


A 


0.35 




3.2 


1.2 


0.27 


0.017 


B 


0.01 




3.3 


1.3 


0.27 


0.017 


C 


0.00 


0.05 




3.74 



Fig. 3 H-R diagram showing the targets for TASK 3 (see Table [JJl. 
Evolutionary tracks correspond to Case 3.1 (bottom-right), Case 3.2 
(middle) and Case 3.3 (top-left). On each track filled circles indicate 
the stages A, B, C. 



Preliminary results of TASK 3 h ave been presented 
by Lebreton et al. (2007); Montalban, T heado and L ebretonI 
(l2007l) ; IChristensen-DalsgaardI (l2008c). Advanced compar- 
isons are presented in this volume fLebreton et al.ll2008bl) . 



5 Specifications for tlie tasks 

In order to reduce the sources of discrepancies in the com- 
parison due to sources not relevant in this context it has been 
necessary to specify the initial parameters and the physics 
to be used in the models. Special care has also been given 
to the exchange of data in order to facilitate the exchange 
of models and evolutionary tracks and their use in detailed 
comparisons. 

Here we briefly describe the common specifications used 
for all TASKS. Further details can be found in the documen- 
tation available at the ESTA web-sit^ 

^ at |http: //www. astro ■ up.pt/corot/compmod/taskl I 
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5.1 Initial parameters 

A list of reference values of some astronomical and physi- 
cal constants have been specified to safeguard consistency in 
the comparison of the output from different codes. We have 
also fixed the mixture of heavy elements to be used in the 
calculations. 

a) Physical constants: We have specified the values of the 
physical constants necess ary for the calculat i on of a stel- 
lar model according to ICohen and Tavloil (Il987l) and 
iLidd d 19941) (see Table [3]l- Furthermore a reference set 
of values is provided for the mass of neutron and atomic 
mass of the elements from 'H to '^O. 

b) Astronomical constants: We have chosen the values to 
consider for the global parameters of the present Sun. 
We have adopted 

/?0 = 6.9599 X 10'" cm, 
for the solar radius (lAllenl Il973h . 

L0 = 3.846 X 10^3 ergs-^ 
for the solar luminosity dWillson et al.Lll98^ and 

M0 = 1.98919 X 10^^ g, 

for the solar mass. Here Mq is obtained from (GMq) = 
1.32712438 x 10^'' cm'^s'^ (|Lidelll994t) and 

G = 6.6716823 x 10"^cm"V's"^- 

The value refers to the radius of the model layer 
where T = T^ff = 5777.54 K. 

c) Initial abundances of the elements and heavy elements 
mixture: All models are calculated with the classical 
GN93 solar mixture of heavy elements dGrevesse and Noelsl 
Il993h . As shown in Table s[T] and |2] the mass fractions of 
hydrogen {X), helium (7) and heavy elements (Z) are 
specified for each model. 

5.2 Input physics 

c) 

In order to be able to make basic comparisons between the 
codes, we chose some reference input physics that have the 
advantage to be currently implemented in stellar evolution 
codes. 

a) Equation of State: We chose to us e the OPAL equation of 
state (EOS) in its 2001 version (iRogers and Navfonovl 
I2OO2I) available at the ORA.L Web sit(fl It consists in a 
series of tables and in an associated interpolation pack- 
age. The tables provide a set of thermodynamic quanti- 
ties, i.e. pressure, internal energy, entropy, specific heat 
Xp = idlnP/dlnp)T, Xi = {dlnP/dlnT)p and the 
adiabatic indices Fi , 12/ (il — 1 ) and (is — 1 ), at given val- 
ues of the temperature, density, X and Z. However, there 

* |http : / /www-phys . llnl . gov/Research/OPAL] 



are differences in the handling of the tables by the differ- 
ent codes. Some use the interpolation package provided 
by the OR\L group, others have developed their own 
interpo lation package. In addition, it ha s been pointed 
out by iBoothroyd and SackmannI ((2003) that there are 
some inconsistencies in the OR\L tables with the re- 
sult that the tabulated values of the adiabatic indexes Fi, 
12/ (l2 — 1 ) and (l3 — 1 ) may differ by several per cent 
from the values that can be recalculated from the tabu- 
lated values of P, Cy, Xp and Xt by means of thermody- 
namic relations. Furthermore, in the course of the com- 
parisons undertaken by the ESTA group it has been shown 
by one of us (I.W. Roxburgh) that the tabulated values of 
the specific heat Cy in OR\L tables are incorrect, and 
also acknowledged by the OPAL team that this param- 
eter should better be obtained from the other quantities 
provided by OPAL. For these reasons, some codes only 
draw P, Xp, Xt and Fi from the OPAL EOS tables and 
apply thermodynamic relations to derive the other ther- 
modynamic quantities. 

Opac ities: We chose to u s e the 1995 OPAL opacity ta- 
bles (^Ig lesias and Rogers'. '1996) complemented at low 
temperatures by the Alexander and Ferguson ( 1994) ta- 
bles, and not to include the conductive opacity. However 
slightly different OPAL tables are actually used by the 
different codes. Pre-calculated OPAL tables have been 
made available by the OPAL group either prior to pub- 
lication or later on the OPAL Web site. They have been 
obtained for a "reduced" GN93 solar mixture, i.e. while 
the GN93 mixture is given for 23 elements, the OPAL 
opacity calculation only considers 19 elements but adds 
the abundances of "lacking" elements to the abundances 
of nearby elements (having close atomic number). Also, 
on the OPAL Web site there is the option to calculate 
tables on-line for any desired 19 elements-mixture and 
to pass or not a smoothing filter to reduce the random 
numerical errors that affect the OPAL opacity computa- 
tion. Different groups used either one of these possibili- 
ties. Moreover, some codes use the interpolation routines 
provided by the OPAL group, others designed their own 
routines. 

Nuclear reaction rates: No specifications have been 
given for the choice of the nuclear network. The basic 
pp chain and CNO cycle reaction networks up to the 
'^O (p,a)'^N reaction have generally been used. Some 
models have been calculated with the option/assumption 
that either one or all the light elements (^He, ^Li, ^Be, 
^H) are at equilibrium while other models have followed 
them explicitly. We chose to compute the nuclear re- 
action rates from the an alytical formulae pr ovided by 
the NACRE compilation ( Angulo et all 1 19991) . Different 
prescriptions have been used for the screening. Most of- 
ten, w eak screening has been assumed under the lSalpeterl 
(119541) formulation. In that case the screening factor is 
written 

/ = exp IaziziJ^] , 



Goals and Tasks of ESTA-CoRoT 



7 



Table 3 Some of the physical constants necessary for the calculation of a stellar model. The units are cgs, except where otherwise noted. 



IjOllZllldllll S CUllSLdllL 


(, — 

K — 


1 ^ROASR V 1 0^1*5 

i.JOUUJO A iVJ 


Atomic mass unit 


m„ = 


1.00Uj4Uz X lU 


Perfect gas constant 


M = 


8.3145111 X 10^ 


Electron mass 


nig = 


9.1093897 X 10"^* 


Electron charge 


e = 


1.602177333 x 10"" C 


Planck's constant 


h = 


6.6260755 x 10"^^ 


Speed of light 


c = 


2.99792458 x lO'" 


Radiation density constant 


a = 


7.5659122 X 10-'5 


Stefan-Boltzmann constant 


a = 


5.67051 X 10-5 


Electron- Volt 


1 eV = 


1.60217733 X 10"" 


Atomic weight of hydrogen 


Ah = 


1.00782500 


Atomic weight of helium 


Aue = 


4.00260330 


lonisation potential for H 


Xh = 


13.595 eV 


1"' ionisation potential for He 


XHe = 


24.580 eV 


2^'^ ionisation potential for He 


XHe+ = 


54.403 eV 



where z\ and zi are the charges of the interactin g nuclei, e) 
Some codes have used the expression (4-221) of lClavtonI 
(Il968h giving A = 1.88 x 10^ and t, = Lz,(l+z,)x, 
where x, is the abundance per mole of element i. Finally, 
in the nuclear reaction network the initial abundance of 
each chemical species is split between its isotopes ac- 
cording to the isotopic ratio of nuclides for which values 
have been given (see the ESTA Web sittH). 
d) Convection and overshooting: W e chose to use the clas - 
sical mixing length tre atment of iBohm-Vitense (fl958h 
under the formulation of lHenyey. Vardya. and BodenheimeJ 
(Il965h taking into account the optical thickness of the 
convective bubble. The value of the mixing length pa- 
rameter has been chosen to be CCmlt = 1 -6. The onset of 
convection is determined according to the Schwarzschild f) 
criterion (Vad— Via(i<0) where Vad and Viad are respec- 
tively the adiabatic and radiative temperature gradient. 
In models with overshooting, the convective core is ex- 
tended on a distance l^y = (Xov x min(//p,/?cc) where Hp 
is the pressure scale height and R^c the radius of the con- 
vective core. We have chosen the value of the overshoot- 
ing parameter to be ttov =0.15. The core is mixed in the 
region corresponding to the convective and overshooting 
region. In the overshooting region the temperature gra- 
dient is taken to be equal to the adiabatic gradient. Note 
that some codes have provided models including over- 
shooting but following different parametrizations (e.g. 
STAROX takes the temperature gradient to be equal to 
the radiative gradient in the overshooting region while 
FRANEC uses a formulation of overshooting depending 
on the mass). 



Atmosphere: Eddington's grey ^(t) law has been pre- 
scribed for the atmosphere calculation: 



r = reff 



3 / 2 



at |http : / /www, astro .up .pt/corot/compmod/taskl | 



where t is the optical depth. The radius of the star is 
taken to be the bolometric radius, i.e. the radius at the 
level where the local temperature equals the effective 
temperature (t = 2/3 for the Eddington's law). Codes 
manage differently the integration of the hydrostatic 
equation in the atmosphere from a specified upper value 
of the optical depth T to the connexion with the envelope. 
Microscopic diffusion: We only considered the diffu- 
sion of helium and heavy elements due to pressure, 
temperature and concentration gradients and we ne- 
glected radiative accelerations that are not yet fully in- 
cluded and tested in the participating codes. We did 
no t impose the formali s m to be used. As reviewed 
by iThoul and MontalMnI (l2007h . two approaches to ob- 
tain the diffusion equation from the Boltzmann equa- 
tion for binary or multiple gas mixtures can be fol- 
low ed: one is based on t he C hapman-Enskog the- 
ory (IChapman and CowlingL ll970L hereafter CC70) and 
the other on th e resolution of the Burgers equations 
(iBurgersl 11969*, hereafter B69). In both methods, ap- 
proximations have to be made to derive the various co- 
efficients entering the diffusion equations, in particu- 
lar the diffusion velocities which are written as a func- 
tion of the collision integrals. In the stellar evolution 
codes which have participated in TASK 3, different treat- 
ments of the diffusion processes have been adopted. 
The TGEC code follows the CC70 approach. The other 
codes follow the B69 approach but differently: (i) the 
A STEC and FRANEC codes use the simplified formalism 
of iMichaud and Proffitd (I1993L hereafter MP93), (ii) the 
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CLES and GARSTEC codes compute the diffusion coef- 
ficients by so lving Burgers' equations accord ing to the 
formalism of iThouI. Bahcall and Loebl (11994) and (iii) 
CESAM provides two options to compute diffusion veloc- 
ities, one based on the MP69 approximation, the other 
on Burge r's formalism, w ith collision integrals derived 
from Pa guette et alj ( Il986l) . Furthermore, the number of 
chemical elements considered in the diffusion process 
differs in the different codes and some codes consider 
that all the elements are completely ionised while others 
calculate the ionisation explicitly. 

5.3 Output format 

Each code had previously designed its own output format. 
But to facilitate the comparisons we have adapted some of 
the output in order to include all information required for 
detailed code-to-code comparisons of models. In most ex- 
ercises we have adopted the file format prescribed by the 
GONG solar model team which is described in the documen- 
tation files in the ADIPLS package and on the ESTA Web sitfi 
This permits a rather easy treatment of the models which can 
be read, plotted and compared by means of the programme 
tools available in the ADIPLS package. In addition, a tool 
MODCONV and associated documentation have been made 
available on the ESTA Web-site by M. Monteiro to allow the 
direct conversion between the different formats available. 

In all cases ESTA participants have been asked to pro- 
vide ASCII files ready for comparisons which contain, for the 
TASK case considered, the global properties of the model 
(name, mass, age, initial composition, luminosity, photo- 
spheric radius, etc.) as well as a wide range of model vari- 
ables at each mesh point (distance r to the centre, mass in- 
side the sphere of radius r, pressure, density, temperature, 
chemical composition, opacity and several other physical 
parameters including quantities of interest for the compu- 
tation of adiabatic oscillations like the Brunt- Vaisala fre- 
quency). 

Modelers have also provided ASCII files giving the vari- 
ation with time of some global or internal properties of their 
evolution sequences (luminosity, effective temperature, ra- 
dius, central hydrogen content, surface helium and heavy- 
element abundances, radius and mass of the convective core 
and depth and mass of the convection envelope, etc). 

5.4 Oscillations 

Stellar oscillations modelers involved in TASK 2 were asked 
to provide adiabatic frequencies in the range [20,2500] /iHz 
and with spherical degrees £=0,1,2 and 3. To obtain the so- 
lution of the equations, the modelers were asked to adopt the 
following specifications: 

a) Mesh: Use the mesh provided by the equilibrium model 
(no re-meshing). 



b) Outer mechanical boundary condition: Set the La- 
grangian perturbation to the pressure to zero {5P = 0). 

c) Physical constants: Use the same as in TASK 1 and 3. 

d) Equations: Use linear adiabatic equations. 

However, different methods and assumptions are made 
to write and to solve numerically the differential equations 
in the participating oscillation codes. In particular: 

a) Set of eigenfunctions: Either the Lagrangian or the Eu- 
lerian perturbation to the pressure {5P or P') have been 
used which affects the form of the equations. 

b) Order of the integration scheme: Either a second-order 
or a fourth-order scheme has been used. 

c) Richardson extrapolation: Some codes using a second- 
order scheme have the possibility to use Richardson 
extrapolation (IShibahashi a nd Osakil Il98lh to decrease 
the truncation error. The combination of a 2"'^ order 
scheme with Richardson extrapolation leads to errors 
which scal e as where ^ is th e number of mesh 
points fe.g. lChristensen-Dalsgaard an d Mullan, 199^). 

d) Integration variable: Either the radius r or the ratio r/P 
have been used as integration variables. 

e) Physical constants: Some codes use a value of the grav- 
itational constant G slightly different than th e one spec- 
ified i n the other TASKS. As pointed out by .Moya et al.l 
(I2008h . using different values of G in the stellar struc- 
ture equations and in the oscillation equations gives rise 
to inconsistencies. 

Th e conseque nces of these choices are discussed in de- 
tails in lMoya et al.. (.2008.) . 



6 Reference grids 

In parallel to the three TASKS, grids of models have been 
especially calculated with the CESAM and CLES codes for 
masses in the range 0.8—10 A/q and chemical compositions 
[Fe/H]=0.0 and — O.IO. These reference grids have been 
used to locate CoRoT potential targets in the H-R diagram 
in the process of target selection (see Figure [Hi and to st udy 
some 5 Scuti candidates for CoRoT dPorettietal-L 120051) . In 
addition, associated oscillation frequencies have been calcu- 
lated for selected models in the grid either with the ADIPLS, 
POSC or L OSC codes. This materia l is described in this 
volume by Montalban et al.' (2008a"); Lebreton a nd Miclieil 
(2008a) : Marques. Monteiro, and Fernandes (200§^ 

In addition, a database of models of j3 Cephei stars and 
their oscillation frequencies called BetaDat has been de- 
signed in the context of ESTA (Thoul, 2008), although the 
models do not use exactly the same input physics and con- 
stants. It is accessible through a Web interfacq^- 



^ at |http: //www, astro .up .pt/corot/ntools/ 



' |http : // astrotheorS ■ astro . ulg . ac . be/ 
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7 Conclusion 

In this paper we briefly presented the different activities un- 
der taken by the evolution and seismic tools activity (ESTA, 
see iMonteiro et all |2006|) organised under th e responsibil- 
ities of the Seismology Working Group (see Mi chel et al.l 
I2OO6) of the CoRoT mission. A description of the ESTA team 
and TASKS aimed at improving stellar evolution codes and 
oscillation codes is presented. This provides the context and 
the reference for the more detailed works published in this 
volume with the descriptions of the codes and the results of 
the comparisons. 

The seismic study of the Sun (helioseismology) had al- 
ready pushed the need to calculate solar models to a new 
level of accuracy due to the necessity to calculate theoretical 
oscillation frequencies that could match the observed solar 
values. As we now approach the same level of precision for 
other stars in asteroseismology the codes used to model the 
evolution of stars of different masses have now to reach the 
same high precision for stellar regimes very different from 
the Sun. As the physics dominating different regimes in the 
H-R diagram are different a strong effort towards this is re- 
quired. 

In the work developed by ESTA we are pursuing that 
goal. The results of this exercise so far have shown that we 
are able to meet the present requirements. But more work is 
planned to address some pending aspects of the physics that 
must be included in the models with sufficient precision to 
be able to reproduce the observed seismic behaviour. This is 
an open project as the arrival of data from CoRoT and other 
future projects will required our tools to improve further the 
precision of the models in order to test the highly accurate 
data made available for asteroseismology. 
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